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Abstract 

A linear defect, viz. an elastic string, diffusing on a planar substrate tra- 
versed by a travelling wave experiences a drag known as Stokes' drift. In the 
limit of an infinitely long string, such a mechanism is shown to be character- 
ized by a sharp threshold that depends on the wave parameters, the string 
damping constant and the substrate temperature. Moreover, the onset of the 
Stokes' drift is signaled by an excess diffusion of the string center of mass, 
while the dispersion of the drifting string around its center of mass may grow 
anomalous. 

PCAS numbers: 05.60.-k, 66.30.Lw, 11.27.+d 
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I. INTRODUCTION 



Particles suspended in a viscous medium traversed by a longitudinal wave f(kx — fit) 
with velocity v = fl/k are dragged along in the x-direction according to a deterministic 
mechanism known as Stokes' drift 0. As a matter of fact, the particles spend slightly more 
time in regions where the force acts parallel to the direction of propagation than in regions 
where it acts in the opposite direction. Suppose |2| that f(kx — fit) represents a symmetric 
square wave with wavelength A = 1-njk and period Tq = 2ir/fl, capable of entraining the 
particles with velocity ±bv (with < b(v) < 1 and the signs ± denoting the orientation 
of the force). Within one cycle, the particle velocity is positive for a longer time interval, 
Th/2(1 — b)v, than it is negative, Tq/2(1 + b)v; hence, the average drift velocity u = b 2 v. 

Recent studies have shown that thermal diffusion of the dragged particles can affect 
markedly the magnitude of the Stokes' drift. A simple change of variables, x — > y = kx — flt, 
makes it possible to map the Stokes' mechanism into the more conventional problem of 
Brownian motion in a tilted washboard potential the dependence of the drift velocity 
u on the viscous constant a and the temperature T of the suspension medium follows 
immediately ||. No matter what the value of a, the Stokes' drift is characterized by an onset 
threshold that grows sharper and sharper as T is lowered toward zero; at finite temperature 
a residual drift occurs in the forbidden parameter region, due to thermal activation {creep). 
Lately, the Stokes's mechanism has been attracting the interest of an increasing number 
of experimentalists in view of its potential applications to single-electron transport in one- 
dimensional channels f|, light storage in quantum wells ||, optical tweezing of colloidal 
particles @, and transport by capillary waves ||, to mention but a few. 

In this article we consider the possibility that a directed linear defect confined to a 
planar substrate can undergo Stokes' drift as an effect of a longitudinal wave travelling 
perpendicularly to its axis. By means of a convenient change of variable, like in Ref. 0, one 
expects to recast the question at hand into an instance of the more general problem of the 



diffusion of an elastic string on a tilted periodic substrate M- [13]. Such a mechanism has 
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been used extensively in the literature to model the diffusion of dislocations in solids || , the 



characteristic curve of Josephson junction arrays [11], the lubrication processes taking place 



between slipping macroscopic surfaces [T3] and the transport of magnetic flux lines in type-II 



superconducting films 1 15 JIT | . Here, as a major difference between point-like [or (0+l)d] and 
linear [or (l+l)d] defects, we show that: (i) In the limit of infinitely long strings, creep effects 
and critically sharp Stokes' thresholds coexhist at any finite temperature; (ii) Such thresholds 
shift toward lower values with increasing the temperature, due to the configurational string 
fluctuations (or phonon dressing); (iii) The dispersion of the drifting string around its center 
of mass grows anomalously with time. 

The paper is organized as follows. In Sec. II, we determine the dependence of the drift 
thresholds on the damping constant and the temperature of the substrate. The onset of the 
Stokes' mechanism is sharp for infinitely long strings, whereas finite-length corrections may 
be computed numerically with no much effort. In Sec. Ill, the phenomenology of the Stokes' 
drift for linear defects is analyzed in detail with particular attention to the excess diffusion 
of the string center of mass in the vicinity of the threshold and the anomalous dispersion of 
the string around its center of mass. In Sec. IV, we compare the effectiveness of the Stokes' 
drift for (0+1) and (1+1) dimensional defects in solids. 

II. DEPINNING THRESHOLDS 

As anticipated in the Introduction, the dynamics of a directed elastic string dragged by a 
longitudinal wave propagating perpendicularly to the string in its glide plane, is intimately 
related to the diffusion of the very same string on a suitable periodic substrate subjected to 
an external static bias. Therefore, we address first the problem of the depinning of a damped 
string from a tilted washboard potential |TI]|]- with the purpose of better characterizing 



the threshold dynamics. 

The perturbed sine-Gordon (SG) equation 



tt -c (j} xx + u sin (f)= -a<p t + F + ((x,t), (1) 



represents the archetypal model of an elastic string in thermal equilibrium on a periodic 
substrate [fO||I5[] . Here, Cq and characterize the string restoring force and the depth 
of the substrate valleys, respectively; correspondingly, Co is the maximum propagation (or 
sound) speed of a signal along the string and uio is the librational (or plasma) frequency of 
the string itself at the bottom of a SG valley. The coupling of the classical field (j){x, t) to 
the heat bath at temperature T is reproduced here by a viscous term — a<pt and a zero-mean 
Gaussian noise source ((x,t). The damping constant a and the noise intensity are related 
through the noise autocorrelation function 

(C(x, t)((x', t')) = 2akT5{t - t')8(x - x'). (2) 

As the relevant SG substrate potential V[<f)] = cJq(1 — cos0) is tilted by the bias term —F(f>, 
the resulting washboard potential retains its multistable structure only if \F\ < F 3 = Uq 
(static depinning threshold ||). 

The elastic string can glide perpendicularly to the V[4>] valleys through three distinct 
mechanisms: 

(a) Kink drift. - A string bridging two adjacent SG valleys bears at least one geometrical 
kink <f> + (x,t), or antikink (f)-(x, t), with <fi±(+oo,t) — <p±(—oo,t) = ±2tt, 

x — x — ut 



(j)±(x,t) = 4arctan 



exp 



dJl — u 2 /cl 



(3) 



Qj 

In the presence of a positive bias F, the solitons <p± are driven apart to =Foo, so that the 
string eventually advances by one step 2ir. SG solitons are known to behave like Brownian 
particles with finite size d = c /uj , mass M = 8/d and average drift velocity ±Up, with 
pOfl up = 2nF/aMo. This mechanism of string transport works for F ^ 0, no matter how 



low the temperature. 

(b) Thermal nucleation. - A string segment lying along one SG valley is unstable against 
nucleation of thermal kink-antikink pairs Jl^j2^|. Thermal fluctuations ((x,t) trigger the 
process by activating a critical nucleus. For F F 3 an unstable nucleus is well reproduced 
by the linear superposition of a kink and an antikink, each experiencing two contrasting 



forces: a restoring force due to the vicinity of its nucleating partner and a repulsive force 
due to either the tilt F itself, or the preexisting (equilibrium) gas of kinks and antikinks 
fl9|j . For F <C F3 the size of the critical nucleus, obtained by balancing the two forces, is 
quite large, so that its energy is about twice the soliton rest energy E = M Cq. At low 
temperature kT <C E , the nucleation time is thus exponentially long, T N oc exp(2E /kT). 
As a result, on establishing opposite, stationary kink-antikink flows, see also mechanisms (a), 
the string creeps in the F direction with average center of mass velocity |17| u t h = 2imoUF, 
proportional to the density of thermal kinks (or antikinks) |19| 



(c) Depinning from the substrate. - Even at low temperature, where thermal creep 
is negligible, the elastic string may unlock from the SG substrate. In the overdamped 
regime, a ^> u , this happens of course for values of F above the static threshold F 3 ; in 
the underdamped regime a <C u)q, the depinning threshold is determined dynamically ||]: at 
T = 0+, stationary string translation is energetically sustainable above the new threshold 

F 2 = 3.36...(— W (5) 

For F > F 2 , the potential energy gained by the string during one jump forward is larger 
than the energy dissipated, in average, by the viscous force that opposes the process. [In 
the noiseless case T = fTTyi2| , of no relevance to the present study, the string dynamics 
would be characterized by hysteretic loops with upper depinning threshold at F 3 and lower 
re-pinning threshold F\ = (4/7r)(a/uJo)Fs.] On increasing a the depinning threshold grows 
monotonically according to the linear law (|5|), for q < Wo, and levels off at the constant 
value F 3 for a ^ 1.2u . 



For tilt values above the depinning threshold F2, an underdamped string slides downhill 
with average velocity uf, both the velocity of the string center of mass and the lateral motion 
of its kinks and antikinks are modulated periodically in time with period As the 
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upper bound of is F/a [§J[13]], the amplitude of the soliton sidewise oscillations is certainly 
smaller than (1/2) [27r/ (F 2 /a)]co, where F 2 is given by Eq. (|5|) and c is the limiting speed 
of any pulse propagating along the string, and therefore smaller than the soliton size d. 
Moreover, for kT <C Eq the mean kink-antikink distance Hq 1 (T) is much larger than d, viz. 
n (T)d <C 1, see Eq. ([J); as a consequence, the spatial configuration of the falling string is 
only negligibly modulated by its passage over the substrate barriers. The same conclusion 
applies in the opposite regime a ^> u . 

Such a periodic effect can be averaged out by adding explicitly a translational coordinate 
X to the field (f>, namely 



(x,t) -> X(t) + <j>(x,t). 



(6) 



On taking the spatial average 



1 r L / 2 

(. . .) x = lim - / (. . .)dx 

L^oo L J-L/2 



(7) 



[not to be mistaken for the stochastic average (. . .) in Eq. (0)] of both sides of Eq. ([[]) along 
a string segment of length L, one obtains 



where we set 



X = -a{X + u th ) - U'{X, T) + F + -=»?(*), 



(8) 



{<f) t (x,t)) x = X + u th , 



(9) 



(10) 



(11) 



(V'[(p)) x = U'(X,T). 



(12) 



In Eqs. (^|) and ( ID ) we have made use of the reflection symmetry dt(j)± = —dt(p T of the 
4>±(x,t) solutions. Note that for kT Eq the <f>± center of mass xq fluctuates with 



velocity variance (u 2 ) = kT/E <C Cg. In Eq. (|j) u t h denotes the creep velocity of the string 
due to the lateral drift of kinks and antikinks, see mechanism (b) above. 

The spatial average of the random force in Eqs. (|l])-(0) is statistically equivalent to the 
Gaussian noise flTT|), where the overall intensity is inverse proportional to the string length 
L, (r](t)) = and 

(7](t)7](t')} = 2akT5(t - t'). (13) 

The effective potential U(X, T) can be computed only under certain low-temperature 
approximations. Of course, it retains all the spatial symmetries of V [</>], namely U(—X, T) = 
U(X,T) and U(X + 2tt,T) = U(X,T). In view of the rigid translation (^j), the identities 
(y^])^ = c<jQ(sin(X + 4>))x — uj 2 [sm X {cos (f)) x + cosX(sin0) x ] and {sm(p) x = 0, lead to 

U(X,T) = V(X)(cos<f ) } x (14) 

with V(X) = w 2 (l — cosX). [Note that the averages (0) (cos(f>) x and (sin0) x are to be taken 
along the kink profile.] The T— dependent configurational factor (cos<f)) x deviates from unity 
because of (i) the local steps (solitons), (ii) the extended excitations (phonons) sustained 
by thermal fluctuations. The first correction is proportional to the (anti)kink density n Q (T) 
and therefore exponentially small; neglecting the soliton contribution corresponds to replac- 
ing (cos (f)) x with (cost/^x; where ip(x,t) represents the equilibrium phonon gas of a string 



sitting at the bottom of one SG valley. The field ift(x,t) can be regarded flq]T9|j as a lin- 
ear superposition of phonon modes ipk(%, t), with dispersion u>\ = Uq + c^k 2 and average 
amplitude set by the equipartition theorem {ipl) x = kT/uj\; hence 

. (2 . kT r+°o dk kT 

At low temperature the distribution of the phonon field is well approximated by a Gaussian 
function P[<f>] = exp(-^ 2 /2(ip 2 ) x ) /J2ir(^ 2 ) x , so that 



(COS ( 



(j,) x = 5ft J / e ^P[iP}diP \ = e~ { ^ /2 . (16) 
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On combining Eqs. ([TJ]) and (|16D we conclude that 

U(X,T) =cj 2 (T)(1-cosX) (17) 

with 

u(T) = u e- kT ' Eo - (18) 

hence, the T dependence of the string depinning threshold (H). 

Let us discuss now the diffusion of the string center of mass described by the reduced 
Eq. (H) or, after the Galileian transformation Y = X + Utht, by 

Y = -aY - u 2 (T) sin(F - u th t) +F+ . (19) 

vL 

The processes X(t) and Y(t) represent a driven Brownian motion on a tilted washboard 
potential ||. In the limit of an infinitely long string the effective temperature of the process 
(pf),(|l9"D T* = T/y/L tends to zero. This implies, for instance, that the depinning threshold 
of an underdamped string is well approximated by Eq. (|5p after u was replaced by u>(T). 
As shown in Risken's textbook |IJ], the depinning mechanism is abrupt; for F = F 2 (T)+ the 
average drift velocity of the string u$ = (X) jumps all of a sudden from its creep value u t h 
up to close to its asymptotic value F 2 (T)/a. In Eqs. (§) and ( (Of) the dependence on the 
string temperature T is restricted to the effective potential amplitude (0) [as a result the 
threshold value F 2 {T) gets lowered] and to the creep velocity u t h{T). Note that even the 
dispersion of the creep velocity tends to vanish in the limit L — ► oo, as the creep mechanism 



through pair nucleation is essentially a shot noise process [pT |. 

Finite-length corrections to the picture outlined above can be computed explicitly by 
applying the numerical recipes of Ref. for a finite effective temperature T*. Depinning 
curves for vanishingly small T* values are displayed in Ref. anticipation of the depinning 
thresholds with the string temperature is hinted at, though not explicitly discussed, in Refs. 

MS 
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III. STOKES' DRIFT 



We go back now to the initial problem of a directed elastic string (say, a lattice disloca- 
tion) confined to a planar substrate and subjected to a longitudinal wave f(<f> — vt) (a stress 
wave) propagating perpendicularly to the string direction. Let the string be coupled to a 
viscous heat bath at temperature T; the corresponding field equation reads 

<Ptt - cl<p xx = -a<p t + f{<f> - vt) + C(ar, t), (20) 

where ((x, t) is defined in Eq. (||) and v = fl/k denotes the propagation speed of a travelling 
wave with wavelength 2ir/k. By means of the Galileian transformation <p — > ip = <fi — vt, Eq. 
(f20p can be rewritten as 

4>tt - c 2 ip xx = -aip t + f{ip) -av + ((x, t). (21) 

For a sinusoidal wave with k = 1, f{ip) = —f^sinip, this equation coincides with Eq. (|]) 
under the correspondence rule 

-> fo, -av, (22) 

and a,T and c unchanged. 

A. Onset thresholds 

The center of mass of the string ip(x,t) diffuses according to Eqs. and (|19D, after 
transformation (p2[) . As long as i/)(x, t) is locked to the bottom of a (creeping) SG valley, see 
Eq. fll9l), the initial string <f>(x,t) drifts with velocity ~ v; upon depinning, the center of 
mass of tp(x,t) travels to the left with negative velocity ~ F/a = —v; correspondingly, 
<p(x, t) gets traversed by the propagating wave without undergoing any appreciable drag, 
were not for the residual creep velocity u t h- 

The relevant Stokes' threshold vg admits two simple limiting approximations; in the 
over damped regime a ^> u 
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v s {a -> 00) = v 3 = — = , (23) 

a a 

in the under damped regime a <C u 

Vs ( a ^0)=v 2 = ^ = 3.36 . . . f(T), (24) 

where f(T) = f exp(—kT/E ). Stokes' drift requires that v < v$, i.e. travelling waves 
must be selected with sufficiently low frequency or short wavelength. At low damping the 
threshold velocity v$ — v 2 is insensitive to a and the drift effect is more pronounced. In the 
limit of infinitely long strings, L —>■ 00, the threshold mechanism grows sharper and sharper, 
as to be expected for the diffusion of a Brownian particle in a tilted washboard potential at 
zero temperature T*. At variance with the (0+1) dimensional case 0,B], the onset of the 
Stokes' drift for an elastic string in thermal equilibrium is thus a critical transition, even at 
finite (but relatively low) temperature. 



B. Excess diffusion 



As proposed in Ref. ||23|| , a dynamic threshold can be characterized in terms of the 
diffusion coefficient 

J= Urn -<*('»' (25) 

of the process at hand. Here, the function D(v) tends to an exponentially small value 
D(0) oc exp(-kT*/f(T) 2 ) for v -> 0, and to the Einstein constant D(oo) = kT*/a for 
v — > 00 B. Both limits vanish for an infinitely long string, T* = T / \f~L — > 0. More notably, 
in the vicinity of the Stokes' threshold v$, D(v ~ v$) is expected to shoot up drammatically; 
such an effect, know as excess diffusion, may be enhanced either by lowering the temperature 
or by weakening the viscous coupling, independently. In the overdamped regime (better said, 
in the Smoluchowski limit a — > 00) and for T* — > 0, the diffusion constant D(v 3 ) decays 
much slower than D(oo), that is |24j 



10 



^3) „ Q 



f 2 (T) 
kT* 



ocLa, (26) 



with Go — 0.90. In the under damped regime, the ratio D(v2)/D(o6) turns out to grow 
even faster, although an appropriate scaling law has not been determined, yet |[23|| . In both 
instances the excess diffusion peak is predicted to be centered very narrowly around the 
relevant threshold velocity vg- 

C. Anomalous dispersion 

Away from the Stokes' threshold vs, the string either drifts along with the travelling 
wave, if v < v$, or creeps forward in the same direction, but with marginal velocity u t h oc 
vexp(— E /kT), if v > v$- The spatial configuration of (f)(x,t), or equivalently of ip(x,t), 
is markedly different in the two opposite dynamical regimes. The creeping string <p(x, t) 
corresponds to the case of the string ip(x, t) falling down a tilted periodic substrate ( pip (v > 



vs)'- As shown in Ref. ||13|| , under stationary periodic boundary conditions the dispersion 

S(t) = ((ij 2 (x,t)) x -^(x,t))l) J (27) 

of an unlocked string segment of finite length L is insensitive to the substrate, S(oo) = 
kTL/Ac^. Accordingly, the spatial dispersion of a creeping string with L — > oo is expected 
to be the same as for a directed elastic string diffusing in a plane, free of periodic modulations, 
either in space or in time, i.e. 

s{t) = k S(±y. (28) 



Co \7ra, 

The counterpart of a drifting string 0(x, t) is a string ip(x, t) locked to the tilted substrate 
(pl|); therefore, the spatial configurations of both strings result in a linear superposition of 
diffusing kinks and antikinks of size d = c / f , separated by an average distance n ~ 1 (T). 
Here the soliton width depends on the wave amplitude fo and on the string elastic constant 
Co, while the soliton height concides with the wavelength 2tt (for k = 1). The question as 
how the dispersion S(t) of an infinitely long string pinned to its SG substrate evolves with 
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time, due to thermal pair nucleation, has been addressed in Ref. p5|: Under the further 



assumption that kink-antikink collisions are destructive, namely for |T7 

v<2^Jv 3 , (29) 

the creep dispersion grows anomalous with fractional exponent larger than 1/2, that is 

S(t) ~ 0.32(n t ^) 2 / 3 . (30) 

Condition (|29|) , clearly satisfied in the overdamped regime, may apply to lower damping 
values, too, provided that v is taken small enough. 

In conclusion, while for v > v$ the string <p(x,t) is almost transparent to the travelling 
wave f(<f> — vt), for v < Vs the string spreads out anomalously with exponent 2/3 and its 
center of mass drifts with velocity close to v. 



IV. CONCLUSIONS 

Longitudinal waves in solids may drag along point-like and linear defects, alike, a mecha- 
nism known as Stokes' drift. The onset of such a transport process, though, exhibits distinct 
properties depending on the dimensionality of the movable defects involved: 

(i) The depinning threshold is very sharp in the case of a long defect line, no matter what the 
substrate temperature. Increasing the temperature may push indeed the Stokes' threshold 
Vs toward lower values, but causes no broadening of the activation region, at variance with 
the depinning of point-like defects. In view of this remark, we conclude that Stokes' drift of 
linear defects is characterized by a clear-cut threshold signature. 

(ii) Due to their spatial extension, linear defects subjected to Stokes' drift undergo a sub- 
stantial deformation in the form of an enhanced transversal dispersion. Moreover, a dragged 
defect line acquires a multi-kink profile, where the geometry of each individual (anti)kink 
component depends on the incoming wave parameters. The ensuing internal dynamics of 
the transported linear defects adds to the dissipative properties of the propagation medium 
overall [p6fl . 
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